Experimental Study of Parametric Autoresonance in Faraday Waves 
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The excitation of large amplitude nonlinear waves is achieved via parametric autoresonance of 
Faraday waves. We experimentally demonstrate that phase locking to low amplitude driving can 
generate persistent high-amplitude growth of nonlinear waves in a dissipative system. The experi- 
ments presented are in excellent agreement with theory. 
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Introduction. When a nonlinear oscillator is resonantly 
driven by small amplitude periodic forcing, the amplitude 
growth is arrested, even at zero dissipation, when nonlin- 
earity comes into play. This is because a frequency mis- 
match develops between the (constant) driving frequency 
and the (amplitude dependent) oscillator frequency 
Persistent amplitude growth can be achieved, by autores- 
onance, when the system nonlinear ly locks to an exter- 
nally varied ("chirped") driving frequency to retain res- 
onant conditions. The precise form of the chirp is unim- 
portant once its sign is correct, and the chirp rate is 
below a critical value. First predicted for harmonic forc- 
ing, autoresonance has found many applications 2]. The 
technique was extended to weakly dissipative oscillators 
1^1 and to nonlinear waves and vortices in non-dissipative 
systems 0|. The theory of parametric autoresonance 
(PAR) was recently developed, first for nonlinear oscilla- 
tors p and later ^] for nonlinear Faraday waves: stand- 
ing gravity waves on a free surface of a fluid which are 
excited parametrically by vertical vibrations. This the- 
ory predicts that a downward chirp of the vibration 
frequency should cause persistent wave growth, which is 
only expected to terminate at large amplitudes, when an 
underlying constant frequency system (CFS), introduced 
below, ceases to exhibit a non-trivial stable flxed point. 

Here we report the flrst experimental veriflcation of 
PAR excitation of a nonlinear wave. Using Faraday 
waves we demonstrate that autoresonance is not hindered 
by moderate dissipation, and the results compare well 
with an extended version of the theory 6] . We show that 
the predicted (negative) frequency chirp indeed drives 
persistent wave growth, via the PAR mechanism, to am- 
plitudes that surpass the theory's region of validity. 

Theory. The theory of PAR excitation of nonlinear 
Faraday waves 6] is based on an amplitude expansion 
that extends earlier constant-frequency treatments 0, 
to the chirped frequency case. Here we summarize the 
main predictions of Ref. ^ and extend the model by (i) 
introducing a more general form of driving acceleration, 
and (ii) taking a more complete account of dissipation. 

Throughout this Letter we consider a rectangular fluid 
cell of length /, width w and depth h in the x, 7/, and z 
direction, respectively. To avoid three-wave interactions 
M we assume that the surface tension corrections are 



small Furthermore, we assume a deep water limit 

h > I and a sufficiently small w so that the fluid motion 
is two-dimensional, depending on z and t. The vertical 
displacement of the vibrating cell is ({t) = (o{t) cos$(t), 
where ^{t) = j^uj{t')dt' is the driving phase, while the 
driving frequency a;(t) and amplitude Co(^) vary slowly on 
the time scale of o;"^. In the weakly nonlinear regime, 
the (time-dependent) scaled acceleration of the cell is 
£(t) = uj'^{t)(o{t)/g <C 1, where g is the gravity ac- 
celeration. As a result, the wave steepness parameter 
krjmax <^ 1, where k = tt/I is the wave number of the 
fundamental mode, and rj is the wave amplitude. As the 
nonlinear frequency shift of standing gravity waves, in the 
deep-water limit, is negative 7], the PAR driving must 
use a negative frequency chirp: duo/dt < 0. The ampli- 
tude dynamics of the fundamental mode are governed, at 
leading nonlinear order, by the equation 6]: 

r)i+ 27771+ y (57)?7?i-31]V)+ ^^'[l+^(^) COS^(t)]77i =0, 

(1) 

where Q = {kg)^^'^ and 7 <C 1^ are the linear wave fre- 
quency and effective linear damping rate, respectively 
(see Ref. 12j for a review of different contributions to 
7). Higher order modes are enslaved to r]i and can be 
calculated once 771 is found. For concreteness, we as- 
sume a linear chirp: u;{t) = cjq — jat^ where /i > is 
constant, and coq = 2Q is the resonant value of the driv- 
ing frequency. We also assume e(t) = eo{l -\- Pt)^ where 
So > and > are constant Now we make an 

ansatz r]i{t) = A{t) cos [Qt + ^(t)], where A and (p are the 
slowly varying amplitude and phase, and use the averag- 
ing method Rescaling time r = {eoQt)/4^ amplitude 



B kAe, 



-1/2 



chirp rate m = 8/i/(£o^)^ and damping 



rate F = 47/ (£0^)7 and denoting (/)(t) 
^ = 4/3/(60^), we obtain: 



liit^/2^(p{t) and 



B = {1^ pr)Bsm{2(l)) -TB , 
^ = (l + ;5r)cos(2^)-5^+mr. 



(2) 



where the dots stand for derivatives with respect to the 
slow time r. The underlying CFS corresponds to m = 
P = 0. Let us start the frequency chirp from the steady 
state obtained for a constant-frequency driving (which is 
the stable flxed point of the underlying CFS). For small 
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m and /3, the PAR wave growth corresponds to the stable 
quasi- fixed point of Eqs. To leading order 



BUt) = [(l+/3r)2-r2 
Mr) = ^arcsin^^ 



1/2 



(3) 



The PAR breaks down if the rescaled chirp rate m ex- 
ceeds a critical value rricr = 0{1) which depends on F |^ 
and p. In any case, the wave growth must terminate at 
large amplitudes, when higher-order corrections to Eqs. 
(QJ and ((2j cause the disappearance of the non-trivial 
stable fixed point in the underlying CFS 0, • 

Alternatively, we can start from a very small initial 
amplitude Bq far from resonance and pass through the 
resonance. In the linear locking stage, we can drop the 
term in Eq. ((2j and obtain 6]: 



BI{t) ^ B^exp 



, ; -77 arcsin(mr) 



(\)^(t) - + - arcsm(mr) - . 

4 2 4^1- (mr)2 



(4) 



where we have put = 0. As long as 5*(r) <C 1, Eqs. (@|) 
are valid on the time interval — 1 < mr < 1 (but not too 
close to vfiT = 1 6]). Remarkably, Eqs. (EJ correspond 
to an unstable (saddle) quasi- fixed point 6], so the sys- 
tem eventually escapes from this point and either enters 
the nonlinear phase-locking regime, described by Eqs. 
with a shifted time, or loses phase- locking. If/when 
B^{r) approaches unity, Eqs. Q become invalid. The 
amplitude B^{r), given by the first of Eqs. reaches 
a maximum at = Vl — F^/m. Its maximum value 



Bq exp 



^— f arccos F 
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decreases with an increase of m. Therefore, at sufficiently 
large m Eqs. and © remain valid over the whole 
interval — 1 < mr < 1. 

Experiment. Our experiments were conducted in a 
transparent cell mounted on a Unholtz-Dickie model 
5PM electro-mechanical shaker made to oscillate in the 
vertical (z) direction. At coq = 54.7±0.16 sec~^ we excite 
an almost pure gravity wave 10] with k = 27r/8 cm~^, 
whose wavelength is twice the cell length of / = 4 cm. The 
cell, of width w = 2 cm, was filled to depth h Q cm 
with hexamethyldisiloxane, and sealed to prevent evapo- 
ration. Hexamethyldisiloxane is a Newtonian fiuid whose 
kinematic viscosity, surface tension and density are, re- 
spectively, 0.65 cSt, 15.6 dyne/cm and 0.76 g/cc. The 
kinematic viscosity was stabilized to within ±1.5% by 
fixing the fiuid temperature to 26.6 ± 0.2 °C. We gener- 
ated acceleration profiles of the form a{t) cos(cc;ot — /it^/2) 
by computer. a{t) was controlled to 1% and measured to 
0.001 ^ resolution by an ADXL103 accelerometer. 



Computer-triggered visualization of the wave profile 
was performed by uniform illumination of the fiuid-air 
interface from behind. The interface's high curvature, 
due to its wetting of the side walls, refracted light away 
from a CCD camera mounted on the opposite side of the 
cell. This resulted in a sharp dark edge depicting the in- 
terface, see Fig. ^ A reference mark on the cell enabled 
measurement of its instantaneous vertical position. Edge 
detection produced a vector of the interface's location 
relative to the cell position as a function of time. The 
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FIG. 1: Images at different phases of (a) (top) an initial wave 
state at cjq and a — 0.064 g and (b) (bottom) the same state 
after being autoresonantly driven at a chirp rate of ji — 0.2 
sec~^ to ^ = (cjo — — 0.08 and a — 0.113 g. 



scaled amplitude of the fundamental, kA, measures the 
steepness of the wave profile. To enable direct compari- 
son to theory, we needed to isolate the fundamental of the 
interface waveform. To this end, we measured the differ- 
ence in wave elevation between two points, chosen to be 
symmetrical about the center of the cell on the x-axis. 
This eliminates, by symmetry, all of the even harmonics 
of the interface elevation. Although the third harmonic 
is not filtered out, the resulting systematic error of A is 
only about 1% at kA = 0.2, and does not exceed 6% for 
kA = 0.6. We could therefore ignore the third harmonic 
while comparing our measurements with the theory in 
the weakly nonlinear regime. The statistical error in kA, 
as estimated from steady state data, is ~ 1%. The in- 
stantaneous phase mismatch (j) between the driving and 
temporal response was extracted using complex demod- 
ulation 13] of two time series: of the measured reference 
mark on the cell and of the wave elevation. 

In Fig. [2^a) we present the measured critical acceler- 
ations ac{5) for the Faraday instability of the fiat sur- 
face as function of the scaled detuning 5 = {coq — ix})/ix)Q. 
The system undergoes a hysteretic transition 8] at an- 
other critical acceleration, ah{S): the lowest acceleration 
at which the nonlinear wave remains stable. 

Figures I2Ib) andI2fc) depict steady-state (CFS) mea- 
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surements of kA and ^, respectively, as a function of 5 for 
a = const. Note that a = sets the maximum attain- 
able detuning. Until a = ah-, kA increases rapidly with 
increasing 5 (decreasing uo). Beyond this point, waves 
will rapidly decay. 



would attain in steady-state for each instantaneous value 
of a{t) and 5{t)^ see Eqs. (0). At larger values of /i, how- 
ever, no phase locking occurs. Here (j) diverges rapidly 
away from the steady- state curves, and subsequently kA 
rapidly decays to zero. 
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FIG. 2: (bottom) (a) Measured ^-dependence of ac (circles) 
and the hysteretic region bounded by an (triangles). Dashed 
lines depict cubic interpolation. Inset: closeup of the vicinity 
of ^ = 0. (top) Steady- state measurements of kA (b) and 
(c) for different values of a = const (legend) versus 6. Here 
r>0.65. Dashed lines in (b) denote the maximum detuning 
attainable for each value of a. Inset in (c) is a closeup of the 
vicinity of 5 = 0, dashed lines depict the predicted values. 

The wave damping rate 7 can be extracted from mea- 
surements of ttc, since F = 1 at the instability onset. 
With 7 in hand, we can directly compare our measure- 
ments to theoretical predictions with no other free pa- 
rameters. The measured phase difference for ^ = agrees 
within 5% with the predicted value. The slow increase in 
(j) with S in Fig. [2|^c) is due to higher order nonlinearities. 

Our first series of measurements used a linear chirp, 
5 = fit/uoQ^ starting from a steady state wave with a small 
but finite amplitude, 0X00 = 00^ and a = > a^- Impor- 
tantly, our choice of the system parameters precluded the 
excitation of other linear modes during a negative chirp. 
As shown in Fig. El^a), ah{S) is almost linear with S for 
0<^<0.12. To maximize the frequency range of the 
excitation, we ramped a linearly in time, a{5) = ao + 
to make a nearly parallel {a = 0.609 g) to ah{S). 

An example of an autoresonant state achieved in a typ- 
ical "ramping" experiment is shown in Fig. ^ where im- 
ages at the initial and final values of 6 are presented. 
Note the substantial increase in the wave amplitude. In 
Fig. ini(a)-(d) we analyze the dynamics leading from Fig. 
^a) to Fig. n^b) by comparing experimental measure- 
ments of kA and with theoretical predictions, obtained 
by integrating Eqs. Q numerically for several chirp rates 
for the same ramp and initial conditions For small 
values of /i, phase-locking occurs where, as a function of 
6, both kA and (j) quickly converge to the values that they 
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FIG. 3: Chirps in 8 initiated from a fixed point. Measured 
(a) and computed from Eq. (b) values of kA versus the 
detuning ^, starting from ^ = and ao = 0.064g. F 0.9. 
Chirp rates of /x = 0.09 (black), 0.15 (blue), 0.18 (red), and 
0.20 (green) sec~^ converge to the fixed point (circles). kA 
for \i — 0.21 and 0.25 sec~^ (red and blue dash-dotted lines, 
respectively) diverge from the steady state and decay. The 
dashed line in (b) depicts the line of quasi-fixed points from 
Eq. 0- (c) and (d): measured and computed values of 
0, respectively, (e) Experimental (diamonds) and computed 
from Eq. (squares) values of the critical chirp rate 
a function of ao. 

As the theoretical curves were obtained with no free 
parameters, the agreement between experiment and the- 
ory is striking. Although the transients are slightly more 
long-lived in the theory, identical convergence/divergence 
of both the phase and amplitude of the waves to/from 
their steady-state values is observed for all of the values 
of /i used. 

An important prediction of the theory is that above 
an ao-dependent critical chirp rate, /icr, phase- locking is 
not possible. For relatively small values of ao, there is 
good quantitative agreement between the measured and 
predicted values of /icr, see Fig. OJe). Furthermore, ficr 
still exists for larger accelerations. Even though ficr in- 
creases with ao, as predicted, the predicted and observed 
values of /icr systematically diverge with increasing ao- 
This divergence is not surprising, as for ao > 0.07 g the 
phase-locking occurs for kA > 0.4, where we would ex- 
pect the weakly nonlinear theory to become inaccurate. 
These results imply that PAR for /i < ficr persists far 
beyond the region of validity of the weakly nonlinear ap- 
proximation. Note that transient stages of <p in Fig. [31 
trace an envelope, corresponding to the universal trajec- 
tory (a saddle point) obtained at /i = /icr 6]. 
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Our second series of measurements involve "passing 
through" the resonance at cjq before any waves are ini- 
tiahy excited. We begin from a flat state at a negative 
detuning Sinu < (i.e. co > cjq), with a^S = 0) < a < 
dci^init)- We then apply a linear chirp S = jit/uo^) while 
keeping a fixed throughout the experiment. As 5 is in- 
creased, we pass through a region of d where a > adS). 

The dynamics of passing through resonance are demon- 
strated in Fig. ^ where experimental measurements of 
kA as function of 5 are shown for several chirp rates. 
As in Fig. IHl the distinction between phase locking at 
low values of /i and phase unlocked states at high values 
of fi is clear. A closer look at Fig. El^a), however, re- 
veals that the precise value of /icr is undetermined. Two 
different runs with the identical chirp rate of fi = 0.27 
sec~^ have qualitatively different behavior: one decays, 
while the other phase locks into PAR. The difference in 
these two runs stems from the dependence of jHcr on the 
wave's initial amplitude Aq. As the waves evolve from 
noise, we do not have experimental control over Aq. In 
Fig. 2fb) we present fits of the initial stages of both of 
the /i = 0.27 sec~^ trajectories shown in Fig. ^a) to Eq. 
Q, where the sole fitting parameter is the value of Aq. 
In both the phase-locked {kAo = 0.00269 ± 0.00002) and 
unlocked {kAo = 0.0009±0. 00004) runs the experimental 
points are indistinguishable from the theoretical predic- 
tions. Thus, a difference of about 10 /im in Aq (~ 0.1% 
of the final, phase-locked amplitude) is sufficient to de- 
termine the wave's long-time dynamics. Similarly, we 
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FIG. 4: Passing through resonance, (a) Measurements of 
kA for a = 0.168 g, -0.046 < ^ < 0.08. Solid (dashed- 
dotted) lines depict phase-locking (unlocking) for the values 
of \i (sec~^) presented in the legend, (b) Comparison of mea- 
sured (dots) and computed from Eq. (lines) values of kA 
in the initial stages of the two runs with \i — 0.27 sec~^. The 
different initial values of kA{) [obtained by fitting Eq. (0], 
resulting from low-level noise, determine the eventual behav- 
ior, (c) Computed values of kA^ using Eq. (0, corresponding 
to the parameters of (a), (d) Comparison of computed (solid 
lines) and measured (dotted lines) values of kA for kA < 0.3. 



obtained Aq for all of the runs presented in Fig. ^a) 
and, using these values, present the computed functions 
kA versus 6 in Fig. El^c), for the parameter values used in 
Fig. Efa). As in Fig. 01 the theory, which uses no other 
free parameters, is a strikingly good description of the 
measurements, especially for values of kA < 0.15, where 
the weak nonlinearity condition kA <C 1 is well satisfied. 
This is demonstrated in Fig. ^d) where a close com- 
parison of theory and experiment is performed for the 
three runs where phase locking failed. In the runs where 
kA < 0.15 throughout the entire experiment, theory and 
experiment are nearly indistinguishable. The agreement 
deteriorates when kA > 0.2, when the system is no longer 
in the weakly nonlinear regime. 

In summary, one can control nonlinear Faraday waves 
by employing PAR. The PAR technique remains opera- 
tional for moderate dissipation and well beyond the weak 
nonlinearity. It would be interesting to extend it to multi- 
mode regimes and to other examples of nonlinear waves. 
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